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Electron-hole generation and recombination rates for intravalley and intervalley phonon scattering 
in Graphene are presented. The transverse and the longitudinal optical phonon modes (-E2 9 -modes) 
near the zone center (T-pomt) contribute to intravalley interband carrier scattering. At the zone 
edge (K(K')-point), only the transverse optical phonon mode (A^-mode) contributes significantly 
to intervalley interband scattering with recombination rates faster than those due to zone center 
phonons. The calculated recombination times range from less than a picosecond to more than 
hundreds of picoseconds and are strong functions of temperature and electron and hole densities. 
The theoretical calculations agree well with experimental measurements of the recombination rates 
of photoexcited carriers in graphene. 
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Graphene is a single two-dimensional atomic layer of 
carbon atoms forming a dense honeycomb crystal lat- 
ticei^^. The high mobility of electrons and holes in 
Graphene has prompted theoretical and experimental in- 
vestigations into Graphene based ultra high speed elec- 
tronic devices such as field-effect transistors, pn-junction 
diodes, and terahertz oscillator a 5 i 6 ' 7 ^. The performance 
of many of these devices depends on the electron-hole 
generation and recombination rates in Graphene. For 
example, the diffusion length of injected minority carri- 
ers in a pn-junction diode is proportional to the square- 
root of the minority carrier recombination time. The 
threshold pumping levels required to achieve population 
inversion, and plasmon gain, in graphene also depend on 
the carrier recombination ratesi Graphene photodetec- 
tors have also been realized^. The efficiency of graphene 
based visible/IR/far-IR detectors would also be critically 
dependent on the carrier recombination tiniest. It is 
therefore important to understand the mechanisms that 
are responsible for electron-hole generation and recombi- 
nation in graphene and the associated time scales. 

Intraband scattering in graphene due to acoustic and 
optical phonons has been extensively studie d 12 ! 13 ' 14 ? 15 ] 16 . 
In this paper we calculate the electron-hole generation 
and recombination rates in graphene due to intravalley 
and intervalley phonon scattering. We find that the to- 
tal recombination times range from less than a picosec- 
ond to more than hundreds of picoseconds and are strong 
functions of the temperature and the electron and hole 
densities. Near the zone center (r-point), only the trans- 
verse and the longitudinal optical phonon modes {Ei g - 
modes) contribute to intravalley interband scattering. At 
the zone edge (K(K')-pomt), only the transverse optical 
phonon mode (A^-mode) contributes significantly to in- 
tervalley interband carrier scattering and the calculated 
recombination rates are found to be faster than those 
due to the zone center i^g-modes. The theoretical re- 
sults compare well (within a factor of unity) with the re- 
cently reported measurements of the recombination times 



II. THEORETICAL MODEL 

In graphene, the electronic hamiltonian for wavevector 
k in the tight binding approximation is^, 
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where t is the energy matrix element between the p- 
orbitals on neighboring carbon A and B atoms and has 
a value of ~3.0 eV, and f(k) — X] m ex P(* &-^m)- Here, 
d m are the position vectors from carbon atom A to its 
three nearest B atoms (FigQJ. Near the K (or K') point 
in the Brillouin zone, the energy dispersion relation for 
the conduction and valence bands is, E s (k) = shv \Ak\. 
s equals +1 and -1 for conduction and valence bands, 
respectively, the velocity v equals V3ta/2 — 10 8 cm/s, 
and Afc is the difference between k and the wavevector of 
the K (or K') point. The corresponding Bloch function 
can be written as a sum over the p-orbitals of all carbon 
atoms, 
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For the A-atoms, Cj(k, s) equals 1, and for the £?-atoms, 

it equals sf*(k)/\f(k)\ = sexp(« </>((£;))). A suspended 
graphene sheet supports phonon modes with both in- 
plane and out-of-plane atomic displacements. The modes 
with out-of-plane atomic displacements do not couple 
with electronic transitions (to linear order in atomic dis- 
placements)^ and will therefore be ignored in this pa- 
per. The displacements ua{<1) and ub(q) of the A and 
B atoms, respectively, corresponding to phonons with 
wavevector q satisfy the dynamical equation 3 , 
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FIG. 1: LEFT: Graphene lattice. RIGHT: First Brillouin 
zone. 
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FIG. 2: Phonon dispersion relation of graphene obtained from 
considering fourth nearest neighbor interactions in the dy- 
namical equation. 



where, following Saito et. al£, the matrix D includes 
fourth nearest neighbor interactions caused by in-plane 
bond-stretching and bond-bending displacements. M is 
the mass of a carbon atom and u>(q) is the phonon fre- 
quency. The phonon dispersion relations found by solv- 
ing the above equation are shown in Fig.[5J Although this 
technique is not as accurate as the ab initio methods for 
obtaining phonon dispersions^ 7 -, it is adequate for calcu- 
lating phonon eigenvectors at the high symmetry points 
and these eigenvectors will be needed in the calculations 
that follow. The LA and TA acoustic phonon modes near 
the T-point can cause intraband carrier transitions via 
deformation potential scattering 18 . Zone center acoustic 
cannot cause interband carrier transitions because energy 
and momentum would not be conserved. Therefore we 
concentrate on the optical phonon modes. In the tight- 
binding approximation, the matrix clement between two 
Bloch states of the perturbing hamiltonian due to atomic 
displacements associated with a phonon of wavevector q 
can be written as, 
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Note that the electron wavevector k is measured form 
the zone center and the relation in Q is valid for 
both intravalley and intervalley scattering processes. Ac- 
cording to the approximation made here, only those 
phonon modes that cause bond-stretching couple with 
electronic transitions. The value of dt/db obtained from 
experiments and density functional calculations is ^45 
eV/nrr>i^. The atomic displacements u^(f) and ftg(r) 
for each phonon mode can be written in the standard 
form, 
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p is the density of graphene and equals ~ 7.6 x 10 7 
kg/m 2 . It remains to calculate the eigenvectors &a{q) an d 
e B (q) for phonon modes near the T and K (K') points. 
These can be obtained from the solution of the dynam- 
ical equation in ([3]) and must be normalized such that, 
e^(g).e_A(g) + e* B {q).e B (q) = 2. Near the T-point, the 
degenerate LO and TO phonon modes correspond to the 
two-dimensional Ei g representation of the point group 
D§h and their eigenvectors can be written as, 
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At the K (if')-points, the group of the wavevector is 
D^h and the phonon mode with the highest energy corre- 
sponds to the one-dimensional A[ representation of this 
group. At the Xi-point 4ir/3a (0, — 1), the eigenvetor is 
found from to be, 
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The two degenerate modes at the K (K' )-points corre- 
spond to the two-dimensional E' representation of the 
group D^h, and their eigenvectors at the -KTi-point are, 
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FIG. 3: Interband intravalley and intervalley electron scatter- 
ing by optical phonons in graphene. 



Finally, the lowest energy phonon mode at the K (K')~ 
points correspond to the one-dimensional A' 2 representa- 
tion of the group and its eigenvector at the i^i-point 
is, 
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(9) 



The eigenvectors for all other K (K 1 )-points can be ob- 
tained either directly from or from the eigenvectors 
at the i^i-point using symmetry arguments based on the 
corresponding representation of the group D^h- Using 
the calculated eigenvectors the matrix elements in (|4]) 
come out to be exactly zero for the E' and the A' 2 phonon 
modes at the zone edge. Therefore, within the tight- 
binding approximation that only bond-stretching motion 
causes electron scattering, the E 1 and the A' 2 phonon 
modes do not cause electron transitions. Away from the 
zone-edge, the matrix elements in Q are non-zero for 
the E' and the A' 2 modes. However, the net contribution 
to electron scattering is small enough compared to the 
other phonon modes that it will be ignored. 



III. RESULTS 

Intravalley interband carrier scattering is due to the 
zone center LO and TO phonons (E 2g -modes) . Assuming 
the frequency of these modes to be dispersion free and 
equal to wro, the recombination rates -RrLO arL d i?rTOi 
and generation rates GrLO and GrLO (units: l/cm 2 -s) 
can be obtained using the matrix elements given in ^ 
and the final results can be written as, 
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Here, f(E—Ef) are the electron distribution functions for 
the conduction and valence bands and n p h(huJro) are the 
phonon occupation numbers. The product E(hujro — E) 
in the integrand comes from the initial and final den- 
sity of states. Intervalley interband scattering is due to 
the zone edge A^-mode. Assuming the frequency of this 
mode to be dispersion free and equal to wkOi the re- 
combination and generation rates, Rko an d Gko (units: 
1 /cm 2 -s) , respectively, can be obtained in a similar fash- 
ion. The results are identical to those given in (| 10[) and 
(fTT|) except that the prefactors are 9/2 instead of 9/4 (i.e. 
twice as large), and wro is replaced by wko- The total 
recombination rate R and the generation rate G due to 
both intravalley and intervalley phonon scattering are, 
R = RrLO + Rfto + Rko and G = GrLO + GrTO + Gko- 
Since wko < w ro ( w ko ~ 1300 cm -1 and wpo ~ 1580 
cm -1 ), and because the integrated squared matrix ele- 
ment for scattering by the zone edge A^-mode is twice as 
large as the one for each of the zone center i? 2ff -modes, 
electron-hole recombination due to intervalley scattering 
is expected to be faster than due to intravalley scattering. 

It is interesting to look at the energy dependence of 
the recombination lifetime r r (E) of an electron in the 
conduction band. Assuming n p h w 0, r r (E) is, 
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The above expression shows that if the valence band has 
enough empty states available then conduction electrons 
near the Dirac point recombine the fastest since they see 
the largest final density of states in the valence band. If 
the valence band has empty states only near the Dirac 
point then conduction electrons with energy near the op- 
tical phonon energy will have the shortest recombination 
lifetime. The average recombination time r r is defined 
as, 



min(n,p) 
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(13) 



Fig-Slplots the recombination and generation rates due 
to intravalley and intervalley optical phonon scattering at 
T=77K and at T=300K are plotted as a function of the 
electron and hole density (assumed to be equal). Equal 
electron and hole densities are generated, for example, 
when graphene is photoexcited as was done in recent ex- 



periments 
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Near thermal equilibrium electron 



and hole densities, recombination and generation rates 
due to intervalley scattering (i?KO or Gko) are ~2 times 
faster at T=300K and ~50 times faster at T=77K com- 
pared to intravalley scattering rates (Rtlo + -RrTO or 
GrLO + GrTo)- This difference is because wko < w ro- 
As would be expected, interband optical phonon scatter- 
ing is less effective at lower temperatures and at smaller 
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FIG. 4: Electron-hole recombination (SOLID) and generation 
(DASHED) rates due to intravalley and intervalley optical 
phonon scattering in graphene at T=77K and at T=300K 
are plotted as a function of the electron and hole density 
(assumed to be equal). The top curve in each pair corre- 
sponds to intervalley scattering (Rko or Gko) and the bot- 
tom curve corresponds to intravalley scattering (-RrLO + ^?rTO 
or Gtlo + GrTo). 
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FIG. 5: Electron-hole recombination time r r due to intraval- 
ley and intervalley optical phonon scattering in graphene is 
plotted as a function of the electron and hole density (assumed 
to be equal) for different temperatures. 



carrier densities. However, even at low temperatures the 
scattering rates can be very fast if the carrier densities 
are large. Fig. [5] shows the average recombination time 
r r due to phonon scattering plotted as a function of the 
electron and hole density (assumed to be equal) for differ- 
ent temperatures. Fig. [5] shows the recombination time 
r r plotted as a function of the temperature for different 
carrier densities. It can be seen that for carrier densities 
larger than mid-10 11 cm -2 the recombination times are 
shorter than 10 ps even at low temperatures. 



FIG. 6: Electron-hole recombination time r r due to intraval- 
ley and intervalley optical phonon scattering in graphene is 
plotted as a function of the temperature for different electron 
and hole densities (assumed to be equal). 
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FIG. 7: Minority carrier (electron) average recombination 
times due to intravalley and intervalley phonon scattering are 
plotted as a function of the minority carrier density for dif- 
ferent temperatures. The majority carrier (hole) density is 
assumed to be 10 12 cm -2 . 



Another important quantity for device applications is 
the minority carrier recombination time in doped (chem- 
ically or electrostatically) graphene. Assuming electrons 
to be the minority carrier, Figs.[7]and [8]plot the minority 
carrier average recombination times due to phonon scat- 
tering as a function of the minority carrier density for 
different temperatures. The majority carrier density is 
assumed to be 10 12 cm" 2 (Fig. [7J) and 10 13 cm" 2 (Fig.©. 
The temperature dependence of the minority carrier life- 
time is seen to depend on the majority carrier density. 
For small majority carrier densities, the minority carri- 
ers near the Dirac point cannot recombine and therefore 
the minority carrier lifetime decreases with temperature. 
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FIG. 8: Minority carrier (electron) average recombination 
times due to intravalley and intervalley phonon scattering are 
plotted as a function of the minority carrier density for dif- 
ferent temperatures. The majority carrier (hole) density is 
assumed to be 10 13 cm -2 . 
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FIG. 9: The calculated net recombination rate, defined as R — 
G, due to interband phonon scattering is plotted as a function 
of the electron and hole density (assumed to be equal) along 
with the experimental data from George et. alfii. T=300K. 



For large majority carrier densities, the minority carriers 
near the Dirac point can recombine very fast and there- 



fore the minority carrier lifetime increases with temper- 
ature. Fig. [5] shows that the minority carrier recombina- 
tion times can be shorter than a picosecond for highly 
doped graphene. For large majority carrier densities, the 
minority carrier lifetime for small minority carrier densi- 
ties and low temperatures follows from (fT2|) and equals 
r r {E = 0) (w 0.34 ps). 

Recently, electron-hole recombination rates in 
graphene were measured using ultrafast optical-pump 
and terahertz-probe spectroscopy^- Fig. [9] shows a 
comparison of the experimental data with the present 
work. Fig. [9] plots the calculated net recombination rate, 
defined as R — G, due to interband phonon scattering 
as a function of the electron and hole density (assumed 
to be equal) at T=300K. Also shown are the best 
fits to the measured time-resolved data from George 
et. alii. The theory compares well (within a factor 
unity) with the experiments. The measured rates are 
between 2.0-2.5 times faster than the calculated rates. 
This difference is likely due to the exclusion of other 
recombination mechanisms in calculating the rates 
plotted in Fig. [5] For example, Coulomb scattering 
(Auger recombination)^ and plasmon emission! can 
also contribute to electron-hole recombination. In 
particular, Auger recombination rates calculated by 
Rana 2 — are of nearly the same magnitude as the phonon 
recombination rates shown in Fig. [9] for carrier densities 
in the 10 11 -10 12 cm~ 2 range. In addition, the presence 
of disorder can also enhance recombination rates, as 
was experimentally observed by George et. al^, by 
providing the additional momentum needed to satisfy 
conservation rules in interband scattering proceses which 
would otherwise be prohibited. 
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